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ABSTRACT 

In most eukaryotes, small RNA-mediated gene si- 
lencing pathways form complex interacting net- 
works. In the ciliate Paramecium tetraurelia, at least 
two RNA interference (RNAi) mechanisms coexist, 
involving distinct but overlapping sets of protein fac- 
tors and producing different types of short interfering 
RNAs (siRNAs). One is specifically triggered by high- 
copy transgenes, and the other by feeding cells with 
double-stranded RNA (dsRNA)-producing bacteria. 
In this study, we designed a forward genetic screen 
for mutants deficient in dsRNA-induced silencing, 
and a powerful method to identify the relevant muta- 
tions by whole-genome sequencing. We present a set 
of 47 mutant alleles for five genes, revealing two pre- 
viously unknown RNAi factors: a novel Paramecium- 
specific protein (Pds1) and a Cid1 -like nucleotidyl 
transferase. Analyses of allelic diversity distinguish 
non-essential and essential genes and suggest that 
the screen is saturated for non-essential, single-copy 
genes. We show that non-essential genes are specif- 
ically involved in dsRNA-induced RNAi while es- 
sential ones are also involved in transgene-induced 
RNAi. One of the latter, the RNA-dependent RNA poly- 
merase RDR2, is further shown to be required for 
all known types of siRNAs, as well as for sexual re- 
production. These results open the way for the dis- 
section of the genetic complexity, interconnection, 
mechanisms and natural functions of RNAi pathways 
in P tetraurelia. 

INTRODUCTION 

Small RNA pathways have greatly diversified in function 
during the evolution of eukaryotes (1,2) but generally rely 



on the production of 21-28 nt short RNAs (sRNAs) from 
single- or double-stranded RNA precursors. Once loaded 
onto Argonaute-containing effector complexes, sRNAs can 
target complementary transcripts and silence homologous 
sequences by diverse mechanisms, such as posttranscrip- 
tional or transcriptional inhibition or even DNA elimi- 
nation (1,3). In RNA interference (RNAi), short interfer- 
ing RNAs (siRNAs) are produced by endonucleases of the 
Dicer family from long double-stranded RNA (dsRNA) 
(4,5). The initial dsRNA triggers may be products of RNA- 
dependent RNA polymerases (RdRPs), enzymes that also 
act at downstream steps in some organisms to amplify the 
silencing response through the synthesis of secondary siR- 
NAs (6,7). 

Different small RNA pathways often coexist in a single 
organism, involving distinct types of sRNAs. In plants, an- 
imals and fungi, miRNAs and endogenous siRNAs regu- 
late the expression of somatic and germline-specific genes, 
and specific endogenous siRNAs control transposons in so- 
matic cells (2,3). Various organisms can also respond to en- 
vironmental dsRNA, for instance some nematode and in- 
sect species that process dsRNA ingested within food bac- 
teria or entering cells by endocytosis (8). Exogenously in- 
duced RNAi pathways were shown in several species to me- 
diate antiviral defense (9-12). Despite the diversity of trig- 
ger molecules and sRNA responses, the different pathways 
involved can be interconnected by sharing some protein fac- 
tors (3). 

In ciliates, a phylum characterized by nuclear dimor- 
phism, a meiosis-specific class of small RNAs (scnRNAs) 
targets the elimination of specific sequences during de- 
velopment of the somatic macronucleus (MAC) from the 
germline micronucleus (MIC) (13-17). As a result, the 
MAC genome is cleared of transposable elements and de- 
rived single-copy internal eliminated sequences (IESs) (18; 
for review see 19). In addition, siRNA pathways are ac- 
tive throughout the life cycle, as shown in Tetrahymena 
thermophila (20,21) and Paramecium tetraurelia (15,22,23). 
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In R tetraurelia, experimental induction of gene silencing 
revealed that different trigger molecules recruit different 
RNAi factors and produce different classes of sRNAs. It 
was first discovered that transformation of the MAC with 
high-copy transgenes lacking the 3' UTR results in aber- 
rant transcripts on both strands (22,24) and silencing of ho- 
mologous endogenous genes (24,25). Transgene-induced si- 
lencing requires the Dicer protein Deri (15), the Piwi pro- 
teins Ptiwil3 and Ptiwil4 (23), and the RdRP Rdr3 (22). 
Although Rdr3 appears to be catalytically inactive, it was 
implicated in the regulation of surface antigen gene expres- 
sion and in the accumulation of endogenous siRNAs of 
unknown function (22). Transgene-induced siRNAs are 5' 
monophosphate and modified at the 3' end, presumably by 
^-O-methylation (22). 

A second pathway, which allows Paramecium to pro- 
cess exogenous dsRNA, was found to involve a distinct 
but overlapping set of RNAi factors. Feeding cells with 
Escherichia coli bacteria engineered to produce dsRNA, 
equivalent to the 'feeding' technique in Caenorhabditis el- 
egans (26,27), induces silencing of endogenous genes ho- 
mologous to the dsRNA (28). In this pathway, Deri (15) 
and the two RdRPs Rdrl and Rdr2 (22) are required 
for accumulation of two types of siRNAs: primary siR- 
NAs, predominantly 23 nt in length and produced from the 
original dsRNA trigger, and secondary siRNAs, presum- 
ably produced from the target mRNA (15). In contrast to 
the primer-independent, 5' triphosphate secondary siRNAs 
produced by direct RdRP activity in C elegans (29-31), R 
tetraurelia dsRNA-induced siRNAs are 5' monophosphate. 
Three Piwi proteins (Ptiwil3, Ptiwil2 and Ptiwil5) appear 
to be involved in this pathway (23). 

The MIC genome of P. tetraurelia has undergone at 
least three whole genome duplications (WGDs), resulting in 
~40,000 genes in the MAC genome, which is responsible for 
all gene expression. 68% of genes still retain their duplicates 
from the most recent WGD (WGD1), and the conservation 
of protein sequences suggests that most duplicate pairs are 
functionally redundant (32). WGDs may explain in part the 
large number of genes encoding core RNAi factors: eight 
Dicer- and Dicer-like genes (including one WGD1 duplicate 
pair), 17 Piwi genes (two unpublished) (six WGD1 pairs) 
and four RdRP genes (no WGD1 duplicate). Only a few 
of these genes have been assigned any function, making R 
tetraurelia a very interesting model to study the genetic and 
mechanistic complexity of small RNA-mediated pathways. 
Due to the lack of an established method to transform the 
MIC genome, however, reverse genetics approaches have 
been restricted to RNAi thus far, providing only limited 
possibilities for the functional analysis of genes involved in 
RNAi. Targeting a gene involved in its own silencing (recur- 
sive RNAi (33)) is a self-defeating process which cannot be 
completed, and is not conclusive when no effect is observed. 

In this study, we describe the outcome of a forward ge- 
netic screen for mutants deficient in dsRNA-induced silenc- 
ing, based on lethal dsRNA 'feeding', which was expected 
to yield loss-of-function alleles for genes that do not have 
functionally redundant paralogs. A combined strategy, in- 
cluding the development of a powerful method to identify 
mutations by whole-genome resequencing, allowed us to 
characterize 71 RNAi-deficient mutants and yielded a total 



of 47 alleles for three known and two previously unknown 
genes. Analyses of allelic diversity revealed a non-essential 
dsRNA processing machinery but provided evidence for an 
essential RNAi function during the life cycle, and suggested 
that the screen is saturated for non-essential, single-copy 
genes. The dsRNA- and transgene-induced RNAi pathways 
are further shown to share essential protein factors. 

MATERIALS AND METHODS 

Paramecium strains, cultivation and genetic analyses 

Mutagenesis and other experiments were carried out with 
wild-type strain 51 of R tetraurelia, mating type E. RNAi 
mutants were back-crossed either with strain 51 or with 
strains carrying the nd7-l mutation (34) as a genetic marker 
(after three rounds of back-cross into strain 51). Unless oth- 
erwise specified, cells were grown at 27° C in wheat grass 
powder (Pines International Co., Lawrence, KS, USA) infu- 
sion medium bacterized with Klebsiella pneumoniae the day 
before use and supplemented with 0.8 |xg/ml (3-sitosterol. 
Genetic analyses of mutants were carried out according to 
standard procedures (35). Fl phenotypes were recorded for 
each of the two karyonides resulting from the two exconju- 
gants of at least two conjugating pairs, and at least 30 F2 
clones resulting from autogamy of Fl heterozygotes were 
studied. 

Random mutagenesis and screening for RNAi-deficient mu- 
tants 

UV mutagenesis (254 nm) was carried out on several inde- 
pendent batches of 200,000 pre-autogamous cells, accord- 
ing to Cohen (36), with a modified dose of 650 J/m 2 . Cells 
were allowed to undergo two divisions before autogamy was 
triggered by starvation to make MIC mutations homozy- 
gous in the MICs and MACs of progeny. Three batches 
in which autogamy reached >90% of the cells, and show- 
ing a post-autogamous survival rate of 75-80%, were se- 
lected for screening. To isolate mutants deficient in dsRNA- 
induced RNAi, post-autogamous cells were first fed with 
three volumes of Klebsiella medium (~l-2 vegetative di- 
visions) and starved again to eliminate fragments of the 
old MAC (monitored by DAPI staining), which would 
otherwise complement mutations in the new MAC. Cells 
were then washed and transferred to three volumes of N- 
ethylmaleimide- sensitive factor (NSF) 'feeding medium' (E. 
coli producing NSF dsRNA, see below) and grown until 
mass lethality (~36 h, four divisions). Surviving cells were 
isolated and grown in standard K. pneumoniae-bacteicized 
medium. Each clone was replicated for storage and fur- 
ther experiments. RNAi deficiency was verified by feeding 
each clone with dsRNAs from the unrelated genes ND7 
and ND169 (34,37) both required for trichocyst discharge, 
a quantitative phenotype that allows to distinguish full and 
partial RNAi deficiencies (see below). 

Induction of RNAi and phenotypic analyses 

Production of dsRNA in E.coli strain HT115DE3, feeding 
to Paramecium cells, transgene-induced silencing of ND169 
and monitoring of trichocyst discharge phenotypes were 



7270 Nucleic Acids Research, 2014, Vol 42, No. 11 




-A 

£ coli 
NSF dsRNA 

I 

isolate surviving 
cells 




F1 




autogamy 



\ 



RNAi/RNAi rnai/rnai 
m/m m/m 




RNAi/RNAi rnai/rnai RNAi/RNAi rnai/rnai 
+H- +/+ m/m m/m 



• • • Oy 



rna/ pool 

Figure 1. Mutagenesis screen and identification of unknown RNAi fac- 
tors. (A) UV-irradiated cell cultures were starved to induce autogamy, a 
self-fertilization process that leads to complete homozygocity of the MIC 
genome. Re-feeding cells with K. pneumoniae medium allowed develop- 
ment of the new MACs from the resulting homozygous zygote within two 
cell divisions, leading to expression of the mutated genes. Transfer of cells 
into NSF dsRNA feeding medium allowed isolation of potential rnai mu- 
tants, as silencing of the essential NSF gene leads to lethality only in wild 
type cells (28,42,43). (B) To identify mutations in unknown RNAi factors, 
silencing deficient mutants were crossed with the wild type, resulting in het- 
erozygous Fl progeny. After autogamy, F2 clones become homozygous for 
the rnai mutant (gray) or wild type (black) allele. Most other (UV-induced) 
mutations (m) are not linked with the target mutation (rnai) and segre- 
gate independently. F2 clones carrying the rnai allele were pooled for MAC 
DNA extraction and whole genome-sequencing. In this data set only the 
rnai mutation is covered by 100% of SNP-containing reads, whereas all 
other mutations are ideally covered by 50% SNP-containing and 50% wild 
type reads. 



carried out as described previously (22,38). TVSF-depleted 
cells were checked for viability after 36 h, the time at which 
typically all cells were dead or dying in a successful silenc- 
ing. To reveal partial loss-of-function phenotypes, ND169 
dsRNA producing E. coli were typically mixed in 1:5 ratio 
with ICL7a dsRNA producing bacteria. 

Plasmid constructs and re-sequencing of known genes 

To induce silencing of TVSFand of ND169 by dsRNA feed- 
ing or by transgene (pTI-) plasmid constructs produced pre- 
viously (22,28) were used. For silencing of PDS1 and CID1- 
like genes, fragments of the coding region were cloned into 
the plasmid L4440 and dsRNA was synthesized in E. coli 
HT115 DE3 (Supplementary Figure SI). Complementa- 
tion of mutants was achieved by microinjection of the en- 



tire open reading frame and flanking sequences, including 
the endogenous promotor, cloned into the pUC18 plas- 
mid vector. For details of all plasmid constructs see Sup- 
plementary Table SI. For re-sequencing of known genes, 
coding sequences including flanking regions were amplified 
by polymerase chain reaction using a proof-reading poly- 
merase (Phusion® High Fidelity, Thermo Fisher Scientific, 
Waltham, MA, USA). Products were purified by phenole 
extraction and polyethylene glycol 8000 precipitation (7,5% 
final) and sequenced with marginal and internal primers 
(oligonucleotide sequences on request). 



Purification of MAC DNA, whole-genome sequencing and 
identification of SNPs 

RNAi mutants carrying an unknown mutation were back- 
crossed to wild type and one karyonide of each exconju- 
gant Fl clone was brought to F2 by autogamy. Thirty F2 
clones were isolated and tested for RNAi deficiency by feed- 
ing of ND169 or NSF dsRNA. Mutant clones showing si- 
lencing deficiency were raised to large cultures. Immedi- 
ately prior to DNA extraction, vegetative mutant cultures 
of both mating types were fused and MAC DNA was ex- 
tracted from this pool (~3 Mio cells) according to the pro- 
tocol used for DNA purification from developing MACs 
(18). Briefly, cells were concentrated, washed, supplemented 
with three volumes of lysis buffer (0.25 M sucrose, 10 mM 
MgCl 2 , 10 mM Tris pH 6.8, 0.2% Nonidet P-40) and lysed 
at 4°C with a Potter-Elvehjem homogenizer, allowing cell- 
but not MAC lysis. MACs were collected and purified from 
bacteria and cell debris by ultracentrifugation through a 2. 1 
M sucrose cushion. Recovered MACs were lysed and DNA 
was extracted as described (39). 5 |xg of RNaseA-treated ge- 
nomic DNA were used for library construction. DNA was 
sequenced by single-end (75 bp lead length) or paired-end 
(50 or 100 bp read length) strategies using Illumina GAII 
and HiSeq next-generation sequencers (for overview of se- 
quence output for each mutant F2 pool, see Supplementary 
Table S2). The reads were mapped to the reference P. tetrau- 
relia strain 51 genome (18,32) using the Burrows-Wheeler 
Alignment tool (BWA) (40) with default parameters and the 
mapping indexed with SAMtools (41). Custom Perl scripts 
and the SAMtools mpileup function were used to identify 
positions where at least 80% of the calls differed from the 
reference genome, at positions with at least 10-X coverage 
and less than 300-X coverage. The list of single nucleotide 
polymorphisms (SNPs) was then curated to identify candi- 
date mutations for experimental validation. 



Small RNA analysis by northern blot 

Total RNA extraction and small RNA northern blots were 
carried out as described (22). As previously, two adjacent 50 
nt sense oligonucleotide probes matching to the dsRNA re- 
gion were used to detect ND169 siRNAs, as most of the pro- 
duced siRNAs are antisense to the target transcript. Simi- 
larly, endogenous cluster22 siRNAs were detected with two 
adjacent 50 nt oligonucleotide probes, complementary to 
the top strand (relative to transcription of the 5' marginal 
open reading frame). 
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Figure 2. Alleles identified from a screen for mutants deficient in dsRNA-inducible RNAi. Unambiguous null alleles (purple, rnai- phenotype as defined 
in the text) were obtained only for RDR1 and PDS1. The location of premature stop codons resulting from frameshifting mutations is shown by gray dots. 
The 5 / TA boundary of the only IES (76 bp) in the RDR2 gene is mutated to AA in the rdr2-\.2A allele, leading to complete retention of the IES in the 
MAC (Supplementary Figure S8B). Hypomorphic alleles are shown in blue (rnai+/- phenotype). Gene accession numbers are available in supplementary 
data (Supplementary Table SI). Co-segregation of the phenotype with the mutation in F2 was verified for alleles labelled with an asterisk (*). Alleles 
complemented by injection of linear DNA encoding the wild-type gene with its natural flanking regions are labeled with two asterisks (**). 



RESULTS 

Lethal dsRNA feeding allows isolation of Mendelian mutants 
fully or partially deficient in RNAi 

To isolate mutants of the dsRNA-induced silencing path- 
way, we first treated wild-type cells with UV light to induce 
random mutations in the MIC genome. Cell populations 
were then allowed to undergo autogamy, a self-fertilization 
sexual process in which only one of the parental MIC al- 
leles is retained, and made homozygous in the MICs and 
MACs of progeny (Figure 1 A). Cultures were then screened 
by feeding them an E. coli strain producing dsRNA ho- 
mologous to NSF, an essential gene involved in exocytosis 
and membrane traffic (42,43). NSF dsRN A feeding rapidly 
kills the wild type (28), so that only mutant cells deficient in 
dsRNA-induced RNAi would be able to survive. One hun- 
dred and fifty surviving cells were isolated, and the RNAi- 
deficient phenotypes were confirmed by dsRNA feeding tar- 
geting different genes (NSF; ND7 and ND169, two single- 
copy, unrelated genes involved in the exocytosis of secre- 



tory granules called trichocysts (34,37); for details, see be- 
low, Figure 7B). After sorting out non-viable, false positive 
or strongly hypomorphic clones, a set of 79 RNAi-deficient 
cell lines was established (detailed outcome of the screen in 
Supplementary Table S3). 

Seventy-one lines were further analyzed. Of these, 64 were 
fully deficient in RNAi, showing normal growth rates in 
NSF feeding medium and wild-type levels of trichocyst dis- 
charge upon ND7 and ND169 feeding. Seven lines were only 
partially compromised in RNAi, suggesting hypomorphic 
mutations. We reasoned that feeding cells with a smaller 
amount of dsRNA would be more sensitive to partial de- 
fects in the RNAi machinery. Consistently, diluting the si- 
lencing trigger (e.g. ND169 dsRNA bacteria in a 1:5 ratio 
with bacteria producing dsRNA from ICL7a, which en- 
codes a non-essential centrin (25)) revealed stronger phe- 
notypes in hypomorphic mutants, comparable to the full 
RNAi deficiency, although this still led to complete silenc- 
ing in wild-type cells. 
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Back-crossing selected lines to the wild type confirmed 
that RNAi phenotypes were due to Mendelian and reces- 
sive loss-of-function mutations: Fl heterozygotes showed 
a wild-type phenotype, and F2 homozygotes obtained by 
autogamy of the Fls showed the expected 1:1 segregation 
of wild-type and mutant phenotypes upon dsRNA feed- 
ing (Supplementary Table S4A). Although these segrega- 
tion analyses would not be sufficient to dissociate closely 
linked mutations, the total number of UV-induced muta- 
tions in the genome was later estimated to be ~ 30 per irra- 
diated cell on average after autogamy (Supplementary Table 
S5A; for mutations, see Supplementary Table S5B), indicat- 
ing that in most cell lines the RNAi deficiency is due to a 
single mutation. 

The mutant screen reveals known and new genes involved in 
dsRNA-induced RNAi 

To identify the mutations in RNAi-deficient lines, we de- 
veloped a whole-genome sequencing procedure to distin- 
guish the mutation causing the phenotype from other, irrel- 
evant UV-induced mutations, after a single round of back- 
cross to the wild type (Figure IB). Briefly, Fl heterozygotes 
were taken through autogamy, and sets of ~30 homozy- 
gous F2 clones were tested for their capacity to respond 
to dsRNA feeding. Deep sequencing of MAC DNA from 
pools of >15 independent RNAi-deficient F2s will iden- 
tify the RNAi mutation as being present in 100% of reads, 
while other, unlinked mutations show up only in 50% of 
reads on average. This strategy was used in combination 
with complementation tests and resequencing of genes pre- 
viously implicated in the pathway to identify candidate mu- 
tations in all 71 lines. 45 of them contained mutations in 
RDR1, RDR2 or DCR1 (Table 1), confirming the role of 
these genes in dsRNA-induced RNAi. nullnull 

The remaining 26 lines were found to carry muta- 
tions in either of two new genes, for which one allele 
was initially identified by the whole-genome sequencing 
procedure. The missense mutation 1.8 (Figure 2; Sup- 
plementary Table S6) is located in the predicted gene 
PTETG9 1000 13001 (ParameciumDB), encoding a protein 
of the nucleotidyl transferase family (also referred to as non- 
canonical polyA/U RNA polymerases (44)) (Figure 3A). It 
is related to the first described non-canonical polyA poly- 
merase, Cidl of S. pombe, as well as to S. pombe's Cidl2 
(Supplementary Figure S2), the first protein among others 
of that family shown to be involved in small RNA medi- 
ated silencing (45). The gene was therefore named CIDL 
Cidl is closely related to Rdn2 of T. thermophila (Supple- 
mentary Figure S2; Figure 6A), which physically interacts 
with the RdRP Rdr 1 in an 'RdRC complex (46,47). CID1 is 
expressed at constant levels throughout the life cycle (Sup- 
plementary Figure S3). 

Mutation 3.18 lies in a 2-kb single-copy gene on MAC 
scaffold 6 (PTETG600032001) (Figure 2; Supplementary 
Table S6), expressed throughout the life cycle (Supplemen- 
tary Figure S3). Blast and conserved domain searches did 
not reveal any significant hit, even in other available cili- 
ate genomes (Tetrahymena thermophila, T. borealis, T. el- 
liotti, T. malaccensis (Oligohymenophora) and Oxytricha 
trifallax (Hypotricha). In contrast, the gene is present in 
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Figure 3. Two newly identified proteins involved in dsRNA-induced 
RNAi. Protein alignments were performed using the MUSCLE v4 soft- 
ware (78). (A) Cidl is a nucleotidyl transferase. Aspartic-acid residues of 
the predicted catalytic triad (*) required for nucleotidyl transferase ac- 
tivity of Arab idop sis thaliana Heso-1 (52) and Tetrahymena thermophila 
Rdnl and Rdn2 (47) are mutated in cidl-1.6 and cidl-3A. Pt, Parame- 
cium tetraurelia; Tt, T. thermophila; Sp, Schizosaccharomyces pombe; Ce, 
Caenorhabditis elegans; At, A. thaliana; Hs, Homo sapiens. (B) Pdsl is con- 
served in Paramecium species. Examples of missense alleles obtained for 
Pt-Pdsl in the RNAi mutant screen are indicated with black arrowheads 
(for substituting amino acids, see Supplementary Table S6). Pt, P. tetraure- 
lia; Ppr, P. primaurelia; Pb, P. biaurelia; Ppe, P. pentaurelia; Ps, P. sexaurelia; 
Pm, P. multimicronucleatum, Pc, P. caudatum. 



all sequenced species of the Paramecium genus, allowing 
identification of conserved amino acids (Figure 3A; Sup- 
plementary Figure S4). Deduced from its specificity for 
dsRNA-induced RNAi (see below), we named this gene 
PDS1 {Paramecium dsRNA-induced RNAi-specific protein 

These candidate mutations were shown to co-segregate 
with the RNAi-deficient phenotype in each of the F2 clones 
included in the genome-sequencing pools, and microinjec- 
tion of the cloned wild-type genes into the MAC of mu- 
tants (cidl -3. 19 and pdsl -3. 17) rescued the RNAi defi- 
ciency (data not shown). In addition, depletion of these pro- 
teins by recursive RNAi experiments (see below, Figure 7D) 
resulted in an RNAi-deficient phenotype. We conclude that 
CID1 and PDS1 are required for dsRNA-induced silencing. 

The RNAi response to dsRNA from food bacteria is not es- 
sential for viability 

Our screen for RNAi-deficient mutants was unlikely to re- 
veal genes with functionally redundant WGD1 ohnologs, 
unless disruption of one copy led to a dosage effect. For 
single-copy genes, we could expect null alleles of non- 
essential genes and hypomorphic alleles of essential genes. 
Unambiguous null alleles were here defined as those con- 
taining nonsense mutations or frameshifts resulting in pre- 
mature stop codons, due to indels or mutations in intron 
splice sites (48), except in cases where the best part of the 
protein is conserved and where there is experimental evi- 
dence for a partial RNAi deficiency. An allele was catego- 
rized as hypomorphic if the RNAi deficient phenotype was 
only partial and/or dsRNA feeding-induced siRNAs were 
reduced, but still detectable on northern blots. Indeed, the 
diversity of alleles obtained for each of the genes hit ap- 
peared to reflect their importance for cellular viability (Fig- 
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Table 1. Mutants obtained in a screen for dsRNA-induced RNAi deficiency 




Mutants 


Alleles 


Function 


Homologs 


RDR1 


41 


26 


RdRP 


Tt, At, Ce, Sp 


RDR2 


2 


2 


RdRP 


Tt, At, Ce, Sp 


DCR1 


2 


2 


RNase III Dicer 


Tt, At, Ce, Sp, Dm, Mm, Hs 


CID1* 


16 


7 


nucleotidyl transferase 


Tt, At, Ce, Sp, Dm, Mm, Hs 


PDS1* 


10 


10 


unknown 


Paramecium 



The first alleles to be discovered for the two new genes (*) CID1 and PDS1 were identified by whole-genome sequencing of pooled F2 clones (see Figure 

i). 
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Figure 4. Abolished siRNA accumulation in RNAi mutants. (A) Northern blot analysis of ND169 siRNAs revealed failure or strong reduction of accu- 
mulation of exogenously induced siRNAs in all mutants, consistent with previous findings from depletion of Deri, Rdrl and Rdr2 by RNAi (15,22,23). 
This suggests a role of these RNAi factors upstream or within siRNA biosynthesis or stabilization. It is of note that silencing efficiencies can vary from 
one dsRNA feeding experiment to another, possibly due to contamination by traces of the standard food bacterium Klebsiella (showing partial Ampicillin 
resistance) which overgrows dsRNA-producing E. coli. This may explain differences in the total level of dsRNA-induced siRNAs in different wild-type 
cultures. The ND169 probe corresponds to a 100 bp region of the dsRNA and is sense oriented, as antisense siRNAs represent the predominant fraction 
of siRNAs in the dsRNA region (22). (B) Detection of endogenous siRNAs from an intergenic region of scaffold 22 (cluster22) revealed that rdr2 mutants 
are unable to accumulate these siRNAs. The probe is complementary to the predominant fraction of siRNAs (top strand) of this region. The lower panel 
shows hybridization to glutamine tRNA as a loading control. 
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Figure 5. Saturation of the screen for non-essential single-copy genes in- 
volved in dsRNA-induced RNAi. Gene sizes of non-essential RNAi fac- 
tors (A) plotted against the number of identified mutant alleles obtained 
(B). Extrapolation by linear regression determined the mutation frequency 
to 1 in 172 bp of coding sequence. 



ure 2): 41 rdrl mutants were found, representing 26 dif- 
ferent alleles. Among these were eight putative null alleles 
(e.g. rdrl-5.28 or rdrl -5.1). Missense alleles were mostly 
non-conservative substitutions in conserved residues, two 
of which in the putative catalytic core region (rdrl-1.4 
D1021N and rdrl-3.16 D1021Y) (49-51) (Supplementary 
Table S6). Similarly, the 10 alleles obtained for PDS1 in- 



cluded null alleles, suggesting that this gene, like RDR1, is 
not required for viability. Although no unambiguous null 
allele was found among the seven cidl alleles, some of the 
mutations changed highly conserved amino acids (cidl- 
3.4 D68N and cidl-1.6 D70N) shown to be required for 
uridylyl transferase activity (47,52). In all putative null mu- 
tants tested, a complete loss of dsRNA-induced siRNAs 
was observed on northern blots (Figure 4A and Supple- 
mentary Figure S6). RNAi-deficient rdrl,pdsl and cidl null 
mutants are fully viable throughout the life cycle (vegeta- 
tive growth and sexual events), and no other phenotypic 
anomaly was observed. Furthermore, the F2 (Supplemen- 
tary Table S4B) and F3 generations of an rdrl-3.l/cidl- 
1.8 double homozygote showed normal vegetative growth 
in standard conditions and at high temperature (34° C). 
We conclude that the capacity to synthesize siRNAs from 
dsRNA ingested with food is not essential in laboratory 
conditions. 

Nevertheless, the dsRNA-induced RNAi pathway ap- 
pears to be highly conserved in the Paramecium genus, since 
tblastn searches identified at least one copy of the RDR1, 
CID1 and PDS1 genes in all sequenced species of the P. 
aurelia complex, as well as in R multimicronucleatum and 
P. caudatum (Supplementary Figure S5). Essential factors 
playing a role in this pathway in P. tetraurelia (see below) 
are also highly conserved (not shown). This suggests that 
the pathway may fulfill an important function required for 
long-term survival of these species in their natural environ- 
ment. 
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Figure 6. Cid2 is also involved in dsRNA-induced silencing. (A) Phylo- 
genetic relationship of the Cidl-like proteins clustering with Pt-Cidl (all 
Cidl-like genes identified in the P. tetraurelia MAC genome and gener- 
ation of phylogenetic trees, see Supplementary Figure S2 and Table S7). 
(B) Double knock down experiment of Pt-Cidl -like genes and ND169 by 
dsRNA feeding. ND169 was used as a reporter for silencing; the gene was 
considered as silenced when cells showed complete deficiency of trichocyst 
discharge (trie-) and inhibition of silencing was measured as proportion 
of cells in the culture showing partial (tric+/— ) or complete (tric+) rever- 
sion of silencing. Silencing was significantly inhibited upon knock down of 
CID1 and CID2 (P-values < 0.03 (*), Mann-Whitney U test, significance 
level 0.05; n = 3; standard deviation is shown). Phenotypes were recorded 
72 h after the first feeding. Bacteria were fed in equal amounts. Triple 
knock down experiments (CID3+CID4+ND169; CID3+CID5+ND169; 
CID4+CID5+ND169) did not show inhibition of ND169 reporter silenc- 
ing (data not shown). Cross-silencing between CID1 and CID2 genes is un- 
likely, as the maximum length of perfect homology is 12 nt, and 18 nt with 
one central mismatch. (C) Northern blot analysis of associated ND169 siR- 
NAs. The ND169 probe corresponds to a 100 bp region of the dsRNA and 
is sense oriented. The lower panel shows hybridization to glutamine tRNA 
as a loading control. 



The screen is saturated for non-essential single-copy genes 

The large numbers of alleles obtained for three non- 
essential, single-copy genes suggest that the screen is sat- 
urated for this category of genes. For RDR1, PDS1 and 
CID1, the number of alleles is roughly proportional to the 
length of the gene; the relationship indicates that, with the 
number of mutant lines analyzed here, one mutation was 
obtained for every 172 bp of coding sequence on average 



(Figure 5). Thus, there may not be any other non-essential, 
single-copy gene involved in the dsRNA-induced RNAi 
pathway. 



Essential genes may be involved in dsRNA-induced RNAi 

In contrast to RDR1, CID1 and PDS1, only two alleles of 
RDR2 were identified. rdr2-3.1 is an amino acid substitu- 
tion (E859K) within a less conserved region of the RdRP 
domain, and rdr2-\. 24 retains an intragenic IES (MIC- 
specific) located downstream of the conserved RdRP do- 
main region, due to a T-to-A mutation in one of the TA 
boundaries that are required for excision during macronu- 
clear development (53) (Supplementary Figure S8B). Both 
alleles are hypomorphic, as indicated by partial RNAi de- 
ficiency and by partial loss of siRNAs upon dsRNA feed- 
ing (Figure 4 A and Supplementary Figure S6). The greater 
reduction in siRNA levels in rdr2-\. 24 was paralleled by a 
weaker silencing of the ND169 target using dilute dsRNA 
trigger (not shown). Similarly, for DCR1 only two hypomor- 
phic missense alleles were identified. Interestingly, the mu- 
tations are located within the N-terminal helicase (dcrl-5.5 
E87K) or helicase C domains {dcr 1-5.21 H622D). Dcrl-5.5 
shows a stronger RNAi deficiency than dcrl-5.21, judging 
from both phenotypes and siRNA levels (Figure 4A and 
Supplementary Figure S6). These findings suggest that the 
helicase and helicase C domains play a role in the biosynthe- 
sis of dsRNA-derived siRNAs, unlike the helicase domain 
of C. elegans' Dcr-1, which is only required for a subset of 
endo-siRNAs (54). The fact that RDR2 and DCR1 are only 
represented in the mutant collection by a small number of 
hypomorphic alleles suggests that they are essential genes 
for which null alleles would lead to lethality. Thus, RNAi 
may fulfill essential functions in P. tetraurelia. 



Another, putatively essential nucleotidyl transferase is in- 
volved in dsRNA-induced silencing 

Whereas no other PDS1 -related gene could be identified 
in the P. tetraurelia genome, 22 Cidl-like non-canonical 
poly(A)-polymerases were found (Supplementary Table 
S7). According to phylogenetic analyses (Supplementary 
Figure S2), five of them, including Cidl, fall in a com- 
mon branch with Rdnl and Rdn2 of T. thermophila (Fig- 
ure 6A). None of them is an ohnolog to another. To an- 
alyze their function in dsRNA-induced RNAi, we used a 
recursive RNAi approach, i.e. we silenced each of them 
by dsRNA feeding, simultaneously with ND169 as a re- 
porter (22). Apart from CID1, targeting CID2 also pre- 
vented silencing of ND169 by dsRNA feeding (Figure 6B). 
This was correlated with a reduced accumulation of ND169 
siRNAs (Figure 6C). It is of note that silencing of a gene 
involved in its own silencing has only limited efficiency, ex- 
plaining the partial impairment of ND169 silencing and 
siRNA accumulation. The results indicate that, like CID1, 
CID2 is involved in dsRNA-induced RNAi. Although it is 
a single-copy gene, our screen did not yield a mutant allele, 
raising the hypothesis that CID2 is essential for viability. 
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Figure 7. dsRNA-induced RNAi factors sharing functions with the transgene-induced silencing pathway. (A)ND169 transgene construct pTI- carrying a 3 f 
truncated version of the ND169 coding sequence driven by a bidirectional constitutive promotor. Green fluorescent protein (GFP) was monitored as control 
for injection and expression (not shown). (B) Ndl69 is involved in the exocytosis of secretory granules (trichocysts) docked at the plasma membrane. Wild- 
type cells expell trichocysts upon treatment with picric acid (tric+). By microinjection of pTI- into the MAC, strongly silenced clones showing complete 
trichocyst discharge defects were obtained (trie-) (not shown), as well as three clones with moderate silencing effect (tric+/-), c2, c5, and c7. Note that 
clone c2 showed up to 18% of cells with tric+ phenotypes in control cultures (ICL7a feeding). (C) Knock down of RNAi factor genes by dsRNA feeding 
in clones with moderate transgene-induced silencing phenotypes. De-repression of ND169 silencing was measured as the percentage of cells showing a 
complete wild type phenotype (tric+). A significant difference to the ICL7a control was detected upon silencing of CID2, RDR2, RDR3 and PTIWI13 
(P-values < 10e-5 (***), one-way ANOVA, significance level 0.005). Phenotypes recorded 120 h after the first feeding are shown. Note that pTI- induced 
silencing efficiencies varied slightly in different cultures of the same injected clone. (D) The effect of RNAi factor depletion on dsRNA-induced silencing 
was verified in parallel as a positive control by double knock down, mixing RNAi factor and ND169 dsRNA-producing bacteria in equal amounts. The 
percentage of tric+ cells, i.e. the degree of inhibition of ND169 silencing, is shown (mean of two replicates for CID1 and PDS1). ICLla control dsRNA 
feeding does neither inhibit dsRNA- nor transgene-induced silencing; PTIWI13 is involved in both pathways, whereas RDR3 is specific for transgene- and 
RDR1 for dsRNA-induced silencing, as previously described (22,23). 



Putatively essential RNAi factors are shared with the 
transgene-induced silencing pathway 

Although transgene-induced siRNAs are biochemically dif- 
ferent from dsRNA-induced siRNAs, Deri and Ptiwil3 are 
shared between these two pathways. Like DCR1, PTIWI13 
is possibly an essential gene, as mutants were not obtained 
in our screen, probably due to inefficient recovery of hy- 
pomorphic mutants. Rdr3, a factor specifically involved in 
transgene-induced silencing, also seems to be essential, as 
its silencing results in a reduced vegetative growth rate (22). 
To test whether the newly identified genes are also involved 
in transgene-induced silencing, we induced silencing of 
ND169 by injecting a transgene construct (pTI-) truncated 
at the 3 'end of the coding sequence (Figure 7A), as previ- 
ously described (22). Candidate genes were then silenced in 
injected clones by dsRNA, and the ND169 silencing phe- 
notype was monitored (recursive RNAi, as described pre- 
viously (22)). To increase the sensitivity of the silencing re- 
sponse to depletion of potentially involved factors, we used 
injected clones that showed only moderate ND169 silencing 
efficiency (tric+/— ) (Figure 7B). Cells showed a strong de- 
crease in silencing efficiency upon CID2 silencing and, as 



expected, upon RDR3 and PTIWI13 silencing (Figure 7C). 
Surprisingly, using this experimental setup we were also able 
to detect a significant reduction of transgene-induced si- 
lencing efficiency when RDR2 was co-silenced, which had 
not been observed previously on strongly silenced clones 
using the same RDR2 silencing construct (22). In con- 
trast, targeting CID1, PDS1 or RDR1 did not result in re- 
duced silencing of the ND169 transgene, although dsRNA- 
induced RNAi was inhibited in control experiments (Fig- 
ure 7D). This suggests that these genes are not involved in 
the transgene-induced RNAi pathway. The involvement of 
the putative essential gene DCR1 in silencing induced by 
the ND169 transgene construct (as well as ND169 dsRNA) 
could also be confirmed (not shown), as previously shown 
for a promotor-less transgene (15). We conclude that RDR1, 
CID1 and PDS1 are specific for the dsRNA-induced path- 
way, while genes hypothesized to be essential according to 
the above criteria (RDR2, CID2, DCR1, PTIWI13) are in- 
volved in both dsRNA- and transgene-induced RNAi path- 
ways. 
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Figure 8. Sharing of RNAi core factors by non-essential exogenously and 
essential endogenously triggered RNAi pathways in P. tetraurelia. In this 
study, mutants were obtained for genes labelled with an asterisk, and genes 
shared between the transgene-induced (dark gray circle) and the dsRNA- 
induced pathway (blue circle) were found to be essential. Brackets con- 
necting two genes symbolize WGD1 duplicate genes. CID2 was identified 
by recursive RNAi in this study. Other genes were identified previously 
(15,22,23). The RDR2 gene is required for the accumulation of endogenous 
siRNAs (purple) of the cluster22 locus, as previously found for RDR3 (22). 
Note that production of dsRNA from this locus is speculative, and that 
other endogenous siRNAs may require different factors. The figure is only 
meant to depict genetic requirements, not the subcellular localization or 
mechanistic roles of the corresponding proteins. 



Rdr2 mutants are compromised in endogenous siRNA accu- 
mulation and show reduced viability after sexual reproduction 

Factors of the transgene-induced RNAi pathway may have 
essential functions in any phase of the life cycle, includ- 
ing vegetative growth, sexual events and development. We 
checked whether mutants are impaired in the accumula- 
tion of a cluster of endogenous siRNAs of unknown func- 
tion (cluster22) which are produced throughout the life 
cycle from an intergenic locus (15). Northern blot analy- 
sis revealed that their accumulation was completely sup- 
pressed in the rdr2-\. 24 mutant (Figure 4B). Signal quan- 
tification (Supplementary Figure S6) showed only a slight, 
if any, reduction for the rdr2-3.7 allele, consistent with its 
strongly hypomorphic nature. Tested mutations in RDR1, 
CID1 and PDS1 had little or no effect. Surprisingly, dcrl 
mutants showed only slight (if any) reduction in cluster22 
siRNA levels as well, as also observed after targeting DCR1 
by dsRNA-induced RNAi (not shown). Similarly, target- 
ing CID2 did not result in any detectable reduction (Figure 
6C), though this may be due to inefficient depletion by re- 
cursive RNAi. Despite the molecular phenotype evidenced 



here, the rdr2 mutants did not differ from the wild type in 
viability or growth rate during vegetative growth. 

In contrast, rdr2 mutants showed a range of defects dur- 
ing sexual reproduction. Survival of homozygotes for the 
stronger allele rdr2-\. 24 was found to be affected by a long 
starvation period (>3 days), but only when this immediately 
followed autogamy, and not when the same starvation was 
experienced by vegetative clones that had already under- 
gone 7-8 divisions (Supplementary Figure S7). New MAC 
development nevertheless appeared to be cytologically nor- 
mal, as judged from DAPI staining, and the reason for this 
post-autogamous defect remains unknown. 

Conjugation was also impaired in both rdr2 mutants. Af- 
ter a cross to the wild type, the Fl clones derived from one 
or both parents were unviable for 35% (18/51) of conju- 
gating pairs and the frequency of proper formation of het- 
erozygous new MACs seemed reduced (Supplementary Fig- 
ure S8A, Table S4A). The most striking effect was that au- 
togamy of the viable, heterozygous Fls most often (9/10; 
P = 0.01) gave rise to F2 progeny in which close to 100% 
of clones were unviable or showed regeneration of the Fl 
MAC (Supplementary Figure S8, Table S4A). This was ob- 
served for Fls deriving from both parents and suggests a 
dominant-negative maternal effect of the heterozygous Fl 
MAC with incomplete penetrance, because the rare cases 
where the F2 progeny survived gave the expected 1:1 seg- 
regation of rdr2 and wild-type F2 homozygotes, and both 
types of clones survived equally well. The reason why post- 
autogamous lethality is much stronger for heterozygotes 
than for rdr2 homozygotes, and the nature of the cellular 
process affected, remain unclear. Testing excision of some 
IESs in post-autogamous, newly developed MACs (rdr2- 
1 .24) did not reveal any DNA rearrangement defect for ma- 
ternally controlled (scnRNA-dependent) or non-maternally 
controlled IESs (data not shown). Nevertheless, our obser- 
vations indicate that RDR2 has essential functions during 
sexual reproduction. A similar phenotype was not observed 
in any other RNAi mutant. 

DISCUSSION 

A whole-genome sequencing method to identify mutations po- 
tentiates a powerful screen 

In this study, we have used a powerful and straightforward 
mutagenesis screen to identify genes required for the RNAi 
response to dsRNA feeding in R tetraurelia: dsRNA ho- 
mologous to the essential gene NSF efficiently kills the wild 
type, and can be fed to large populations of cells. Muta- 
tional studies have long benefited from the easy genetics 
and short generation time of this organism (55-59), but 
the identification of mutations relied on functional comple- 
mentation and sib-selection sorting of an indexed library 
(60,61), a time-consuming method ill-suited for large mu- 
tant collections. We have designed an efficient procedure to 
identify phenotype-causing mutations after a single round 
of back-cross to the wild type. Other, unlinked mutations 
segregate independently among F2 homozygotes, so that 
whole-genome sequencing of pools of phenotypically mu- 
tant and wild-type F2s will show the RNAi mutation in 
100% and 0% of reads, respectively, while most irrelevant 
mutations should be present in 50% of reads from both 
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pools. In practice, sequencing of the wild-type pool was 
found to be unnecessary; 15-X coverage of pools of ~20 
mutant clones, with 75-bp single reads, was often sufficient 
to narrow down the list of candidate mutations to a single 
one, although 100-bp paired-end sequencing proved more 
efficient. Distinguishing closely linked mutations would re- 
quire more F2s and a higher coverage, but the standard mu- 
tagenesis conditions used (~20% lethality in homozygous 
progeny of irradiated cells) were shown to result in an aver- 
age of ~30 mutations in the 72-Mb MAC genome of viable 
progeny, making cases of tight linkage infrequent. Similar 
approaches have been used in C. elegans and A. thaliana 
(62,63), but here the possibility to generate large numbers of 
entirely homozygous recombinant F2s by autogamy makes 
the method faster and simpler, dispensing with the need for 
extensive back-crossing or pre-mapping. 

Mutations causing full or partial RNAi deficiency identify 
known and new genes 

Mutations were identified in 71 RNAi-deficient cell lines, 
yielding a total of 47 different alleles for only five genes 
(Figure 8). Three of them (RDR1, RDR2 and DCR1) 
were previously implicated in dsRNA-induced RNAi on 
the basis of recursive RNAi experiments. The two new 
genes included the nucleotidyl transferase CID1 and the 
Paramecium-spQciiic gene PDS1. Among the 47 alleles, 14 
were isolated more than once; in three cases the same alle- 
les were recovered from different mutagenized cultures and 
thus originated from different mutation events. This sug- 
gests that our analysis provides a rather exhaustive view of 
the genes and alleles that can be identified by the experimen- 
tal scheme used. 

The latter allowed recovery of some hypomorphic alleles 
(seven out of 47) causing only partial RNAi deficiencies, as 
confirmed by functional tests using dilute dsRNA trigger. 
The rather low frequency of such alleles may be due to the 
~7 vegetative divisions allowed for fully resistant mutants 
in NSF dsRNA feeding medium, before cell isolation: al- 
though hypomorphic mutants were able to survive that step, 
their growth rate was probably greatly reduced. The exper- 
imental scheme would thus have to be modified to yield a 
higher frequency of hypomorphic alleles, the only ones that 
can be obtained for essential genes. 

The dsRNA- and transgene-induced RNAi pathways share es- 
sential genes 

Analysis of the diversity of alleles obtained for each of these 
five genes showed that they fall in two groups. For RDR1, 
CID1 and PDS1, the many alleles recovered included un- 
ambiguous null ones such as nonsense mutations and/ or 
non-conservative changes of conserved residues, which all 
abolished accumulation of siRNAs upon dsRNA feeding. 
Although these mutants were completely deficient in the si- 
lencing response to dsRNA, they did not show any other 
phenotype at any stage of the life cycle, indicating that the 
capacity to synthesize siRNAs from dsRNA ingested as 
food is not required for short-term viability in laboratory 
conditions. In contrast, only two alleles were obtained for 
each of DCR1 and RDR2 and these were hypomorphic, 



only partially impairing the silencing response and siRNA 
accumulation. This strongly suggests that these are essential 
genes in which null mutations would be lethal at some stage 
of the life cycle. In particular, these results establish a clear 
distinction between RDR1 and RDR2, two RdRPs which 
could not be distinguished in previous recursive RNAi tests 
applied during the vegetative phase (22). 

Strikingly, these two groups of genes also differ in the 
specificity of their involvement in RNAi pathways (Figure 
8). We confirmed that DCR1 is also required for transgene- 
induced silencing and showed that the same is true of 
RDR2, although this had not been seen previously. Here, 
recursive RNAi tests (by dsRNA feeding) were made more 
sensitive by using a transformed clone showing only moder- 
ate silencing of the ND169 reporter (as proposed in a model 
for recursive RNAi dynamics (33)). In the same conditions, 
recursive silencing of RDR1, CID1 and PDS1 did not af- 
fect transgene-induced silencing. Thus, the only essential 
genes uncovered by our screen are those also involved in 
transgene-induced RNAi. Two more genes likely belong to 
this group. By testing close paralogs of CID1, we found that 
CID2 is involved in both pathways, as previously shown for 
PTIWI13. The fact that no mutation was recovered in these 
genes is consistent with the idea that they are also essential. 

A third, presumably non-catalytic RdRP gene previously 
implicated in transgene-induced silencing, RDR3, is further 
required for normal vegetative growth rate and accumula- 
tion of cluster22 endogenous siRNAs. The latter differ from 
transgene-induced siRNAs in that they lack a modification 
of the terminal ribose, and interestingly, were also absent 
or reduced in the hypomorphic rdr2 mutants. RDR2 thus 
appears to be required for all known types of siRNAs in 
Paramecium (Figure 8). However, the wild-type growth rate 
of the mutants suggests that its putative essential function 
may not be needed during the vegetative phase. 

In contrast to rdr2, no effect was observed in dcrl mu- 
tants or CID2 knock down cells on cluster22 siRNAs. Al- 
though we cannot completely exclude a false-negative result 
due to the hypomorphic nature of the dcrl alleles and the 
inefficient knock down of CID2, this may indicate that these 
endogenous siRNAs are not dependent on all components 
of the transgene-induced pathway and are thus produced 
by a different pathway (Figure 8). This is supported by the 
observation that they do not carry a 3 f modification. 

RDR2 is involved in sexual reproduction 

The crosses of rdr2 mutants to the wild type revealed de- 
fects in sexual reproduction. Fls showed a high lethality 
rate and occasional failure of the formation of heterozy- 
gous MACs. However, it is unclear whether this was sig- 
nificant in Fls deriving from the wild-type parent (which 
would point to defects in pre-zygotic stages such as karyo- 
gamy). More striking was the frequent unviability of most 
or all F2 clones produced by autogamy of some — though 
not all — viable heterozygotes, regardless of F2 genotype or 
Fl parental origin. The stochastic nature of the effect raises 
the possibility that it is due to defects in maintenance of the 
germline in mutant homozygotes, such as defects in chro- 
mosome segregation during MIC mitosis, which would have 
no impact on vegetative growth of mutants, but would lead 
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to aneuploidy in the Fl and therefore lethality in the F2. 
Such defects have been observed for instance in RNAi mu- 
tants in S. pombe, including the RdRP Rdpl (64,65) and 
the nucleotidyl transferase Cidl2 (66), which are both essen- 
tial for pericentromeric heterochromatin formation (45,67). 
Stochastic loss of MIC chromosomes during vegetative 
growth would make each mutant clone potentially unique 
in its ability to foster viable F2 progeny. 



The genetic complexity of the dsRNA-induced RNAi pathway 

Previous studies have probed the genetic complexity of 
RNAi pathways (68-72). In C. elegans, a genome-wide 
RNAi screen for genes involved in silencing induced by 
RNA hairpin expression identified 90 genes, 54 of which 
are essential for viability (69). In our study, the screen 
seems to be saturated only for the sub-category of non- 
essential, single-copy genes, all three of which are specifi- 
cally involved in dsRNA-induced RNAi. Additional non- 
essential genes are likely to have been missed because they 
retain their WGD1 duplicates, which very often have re- 
dundant functions. Indeed, 68% of all P. tetraurelia genes 
still have their WGD1 duplicate (32). Such may be the case 
of PTIWI12 1 PTIWI1 5 previously shown to be involved 
in dsRNA-induced silencing but for which co-silencing of 
both duplicates was required to detect an effect (23). Other 
genes may have been missed because they are essential, since 
our scheme yielded only a low frequency of hypomorphic 
alleles. Essential genes may be shared with the transgene- 
induced pathway; in addition, processing of dsRNA from 
food bacteria could be coupled to other cellular processes 
such as uptake of nutrients, as shown in D. melanogaster 
S2 cells and C. elegans (73,74). Thus, the total number 
of genes contributing to dsRNA-induced RNAi might be 
much higher than the nine known so far. 



Environmental RNAi in Paramecium 

Environmental RNAi, i.e. the processing of exogenous 
dsRNA into functional siRNAs, has been observed in a 
number of organisms (8), but has not yet been reported 
to be essential for cellular viability. It is also dispensable 
for short-term viability of P. tetraurelia in laboratory con- 
ditions. However, the genes involved are highly conserved 
in the Paramecium genus, suggesting they may have impor- 
tant functions for the long-term survival of these species 
in their natural environment. Exogenous RNAi may serve 
as an antiviral defense, a major function of the pathway in 
D. melanogaster, C. elegans and other species (9-12), even 
though no Paramecium virus has been identified so far. If 
most of the genes in the pathway are conserved for the pro- 
cessing of viral RNA, only a few more may be required for 
the uptake of dsRNA from phagosomes, and these are not 
necessarily conserved. Such a situation is seen with some 
Caenorhabditis species that lack the dsRNA uptake mecha- 
nism but respond to exogenous dsRNA when the C. elegans 
dsRNA transporter SID-2 is expressed (75). The identifica- 
tion of novel genes and the availability of mutants will help 
investigations into the mechanism and natural functions of 
dsRNA uptake. 



Biosynthesis of dsRNA triggers and siRNAs 

Although it is still unclear how RdRPs are specifically re- 
cruited to their target template during biogenesis or main- 
tenance of silencing triggers in vivo, recent studies in S. 
pombe, C. elegans and T. thermophila point to the forma- 
tion of RdRP complexes (RdRC) as an important step 
(45,47,76,77). In the ciliate T. thermophila, a single RdRP 
(Rdrl) interacts in a mutually exclusive manner with ei- 
ther of two nucleotidyl transferases, forming an essential 
RdRC with Rdnl, or a non-essential one with Rdn2. The 
latter contains an associated factor, either Rdfl or Rdf2 
(47,77). The role of these specific RdRCs in vivo is still 
unclear (21,47). In P. tetraurelia, it is tempting to specu- 
late that Rdrl and Cidl form a non-essential RdRC which 
specifically acts in dsRNA-induced silencing, while Rdr2 
and Cid2 form an essential RdRC acting in both path- 
ways (Rdfl and Rdf2 are not conserved in the P. tetrau- 
relia genome). The reason why two RdRCs should be re- 
quired for dsRNA-induced silencing, however, will remain 
puzzling until their precise functions are identified. They 
may have distinct roles in the production of primary and 
secondary siRNAs (15,22). The present work opens the way 
for an analysis of siRNAs and their precursors in each mu- 
tant, which is likely to provide insight into the processing 
mechanisms. 
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